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ABSTRACT 

We explore the self-enrichment hypothesis for globular cluster formation with respect 
to the star formation aspect. Following this scenario, the massive stars of a first stellar 
generation chemically enrich the globular progenitor cloud up to Galactic halo metal- 
licities and sweep it into an expanding spherical shell of gas. This paper investigates 
the ability of this swept proto-globular cloud to become gravitationally unstable and, 
therefore, to seed the formation of second generation stars which may later on form 
a globular cluster. We use a simple model based on a linear perturbation theory for 
transverse motions in a shell of gas to demonstrate that the pressures by which the 
progenitor clouds are bound and the supernova numbers required to achieve Galactic 
halo metallicities support the successful development of the shell transverse collapse. 
Interestingly, the two parameters controling the metallicity achieved through self- 
enrichment, namely the number of supernovae and the external pressure, also rule the 
surface density of the shell and thus its ability to undergo a transverse collapse. Such a 
supernova-induced origin for the globular cluster stars opens therefore the way to the 
understanding of the halo metallicity distributions. This model is also able to explain 
the lower limit of the halo globular cluster metallicity. 

Key words: globular clusters: general ~ Galaxy: halo - supernova remnants - stars: 
formation. 



1 INTRODUCTION 

Globular clusters (GC) are dense, massive and round-shaped 
groups of stars present in the vast majority of galaxies. In 
our Galaxy, the halo GCs were among the very first bound 
structures to form and their study provides therefore valu- 
able information about the early Galactic evolution. Their 
formation is an exciting but yet unsolved problem. For in- 
stance, it is still an open question whether Galactic halo 
GCs formed out of gas already chemically enriched {pre- 
enrichment models, e.g. Harris & Pudritz 1994) or whether 
they produced their own heavy elements through an ear- 
lier generation of stars within the GC progenitor itself {self- 
enrichment models). In the second class of models, the is- 
sue of their formation is directly related to the origin of 
their metal content (Cayrel 1986; Brown, Burkert & Truran 
1995; Parmentier et al. 1999). Such a feature makes the self- 
enrichment scenario especially appealing if we hypothesize 
that the GC progenitor clouds are made of primordial (i.e. 
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metal-free) gas, a reasonable assumption for the Old Halo 
GCs, i.e. the population of old and coeval (Rosenberg et 
al. 1999) halo GCs. 

Regarding the origin of the proto-globular cluster clouds 
(PGCC), Fall & Rees (1985) suggested that they formed out 
of the collapsing protoGalaxy, as cold and dense clouds in 
pressure equilibrium with a hot and diffuse background. In 
the frame of this theory, the GC progenitor clouds are ther- 
mally supported (i.e. no additional support against gravi- 
tation due to magnetic fields and turbulence) and made of 
primordial gas. In order to explain the metallicities of halo 
GCs, Parmentier et al. (1999, hereafter Paper I) expanded 
the Fall & Rees (1985) model for GC formation by the self- 
enrichment picture. According to this one, a first stellar gen- 
eration forms in the central regions of each PGCC. When 
the massive stars explode as Type II supernovae (SNell), 
they chemically enrich the surrounding gas and sweep the 
cloud, turning it in an expanding shell of gas in which the 
formation of a second, chemically enriched, stellar genera- 
tion may be triggered. These second generation stars form 
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the proto-globular cluster. 

The sweeping and compression of the interstellar 
medium by massive star explosions is not the sole mech- 
anism invoked to account for star formation during the ear- 
liest stages of the Galactic evolution. Some other models 
assume a different origin for the trigger. For instance, fol- 
lowing Vietri & Pesce (1995), the high pressure confining the 
GC gaseous progenitor leads to the propagation of a strong 
shock inwards the cloud, stimulating thereby the formation 
of new stars. On another side, Murray & Lin (1992) and 
Dinge (1997) suggested that the propagation of shock waves 
is promoted by cloud-cloud collisions. In a similar way, star 
forming clouds may have coalesced into larger units until 
they reach a density and/or accumulated mass high enough 
to enable the formation of bound globular clusters (Larson 
1988, Smith 1999). Obviously, several processes are able to 
collect ambient gas into dense layers/clouds where star for- 
mation can thereafter take place. Several of them may have 
been at work at the same time. Our interest being in un- 
derstanding the origin of the metal content of GCs, in what 
follows, we investigate the hypothesis of star formation in 
gas layers swept by the explosions of PopIII massive stars 
located at the center of the PGCCs. 

The debate of how PopIII stars looked like has been rag- 
ing for many years. Numerous studies have addressed the 
issue of the collapse and fragmentation of primordial gas 
clouds in order to estimate the masses of PopIII stars. The 
primordial gas being deficient in heavy elements (i.e., the 
most efficient coolants in present-day star forming clouds), 
all of them have emphasized the importance of cooling 
by H2 molecules. Inspite of this, the achieved conclusions 
do not necessarily converge. Recent numerical simulations 
(Abel, Bryan & Norman 2002) suggest that metal-free stars 
form in isolation and are massive (30 < M < 100 M©) ob- 
jects. From their own simulations, Bromm, Coppi & Lars- 
son (1999) quoted a characteristic mass even larger than 
100 Mq. These results are in marked contrast with the early 
study performed by Palla, Salpcter & Stabler (1983) follow- 
ing which primordial gas clouds should be capable of frag- 
menting into low-mass stars (i.e., down to 0.1 Mq). Naka- 
mura & Umemura (1999) reached a somewhat intermediate 
conclusion. According to them, the mass range of the first 
stars is similar to its present value (i.e. no star with M > 
100 M0) but nevertheless excludes low-mass long-lived stars, 
that is, the lowest mass star allowed to form in a metal-free 
medium is a ~ 3 M© star. In a more recent model, Nakamura 
& Umemura (2001) predict a bimodal initial mass function 
for PopIII stars. In fact, the initial mass function of metal- 
free stars would show two distinct peaks at ~ 1 Mq and ~ 
100 Mq, that is, it would include intermediate mass (1 < M 
< 10 Mq) stars, massive (10 < M < 100 Mq) stars as well as 
very massive (M > 100 M©) ones. As quoted by Christlieb 
et al. (2003), the recent discovery of HE 0107-5240, the most 
metal-poor ([Fe/II] =— 5.3) star ever discovered in the Galac- 
tic halo, could be a challenge to these models precluding the 
formation of stars with mass low enough for their life dura- 
tion to exceed a Hubble time. 

While there is currently a wide consensus that present 
star formation operates predominantly in a clustered mode 
(Lada, Strom & Myers 1993), the question whether clus- 



ters of metal-free stars managed to form in the first (proto-) 
galaxies is still open. In this paper, we assume that at least 
some of the primordial star formation sites produced stellar 
clusters whose most massive stars (10 < M < 50 Mq) end 
their life as canonical SNell. Therefore, our study does not 
include the possibility of a prompt enrichment of the pri- 
mordial interstellar medium by a population of very mas- 
sive (i.e., > 100 Mq) stars as suggested by, e.g., Wasserburg 
& Qian (2000). The chemical enrichment provided by iso- 
lated (very) massive objects will be qualitatively discussed 
in Sect. 3 where we will see that they may be appropriate 
to account for the formation of very metal-poor stars, i.e., 
stars more metal-poor than the most metal-deficient halo 
GCs ([Fc/H] < -2.5). 

Supernovae having long been thought to disrupt the 
cloud of gas out of which they have formed, Parmentier et 
al. (1999) studied the ability of pressure-truncated clouds to 
retain SNII ejecta. This ability will rule the metal content 
of the proto-cluster. In fact, in this class of models, the final 
metallicity is determined by the number of supernovae {N, 
which, assuming a given initial mass function and given SNII 
yields, determines the amount of metals dispersed within the 
PGGC) and the background pressure {Ph, which determines 
the mass of the cloud, that is, the mass of primordial gas 
to be chemically enriched). Comparing the gravitational en- 
ergy of the PGCC with the kinetic energy of the supershell 
resulting from the SNII explosions, Parmentier et al.(1999) 
showed that the gaseous GC progenitors can sustain up to 
~ 200 supernovae (the disruption criterion, Eq. 14, Paper I) . 
This is a number high enough for the PGCCs to achieve halo 
metallicities. Furthermore, for a given number of exploding 
massive stars, a self-enrichment episode in pressure-bound 
clouds lead to a metallicity gradient throughout the resulting 
system of GCs and to a correlation between the mass and 
the achieved metallicity of the gaseous progenitors in the 
sense that the least massive clouds are the most metal-rich. 
Such a trend emerges because if the bound pressure is higher, 
the mass of the pressure-truncated cloud will be lower, and 
its ability to retain supernova ejecta will be greater. These 
trends, i.e. a metallicity gradient and a mass-metallicity re- 
lation, are indeed statistically present in the Old Halo, that 
is, the halo from which the presumably accreted component 
has been removed (Parmentier et al. 2000, Parmentier & 
Gilmore 2001). 

After having analysed to which extent the Galactic halo 
GC data fit the correlations induced by the self-enrichment 
process, the next step is to wonder whether there are some 
stars forming out of the chemically enriched supershell, i.e. 
whether there is a second stellar generation tracing these 
correlations. The idea of star formation induced by super- 
nova explosions dates back at least to Opik (1953). Nu- 
merous cases of distant star forming loops, connected with 
shells resulting from supernova explosions and located in the 
Galactic disc as well as in dwarf galaxies, are detailed in the 
literature (e.g. Comeron & Torra 1994, Walter et al. 1998, 
Efremov & Elmegreen 1998, Rubio et al. 1998). In this pa- 
per, we address the specific case of supersliells made of swept 
PGCCs and within which the formation of the stars of future 
GCs may be triggered. In a first step, we limit our study to 
the propagation of the shell throughout the hot protogalac- 
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tic background in which the PGCCs arc initially embedded, 
this part of the shell propagation being much longer than 
the propagation throughout the cloud (see Fig. 1). 

The outline of the paper is as follows. In Sect. 2, we 
solve the perturbed equations of continuity and motion for 
transverse flows within the shell (i.e. the swept cloud) in or- 
der to identify the conditions supporting a successful shell 
transverse collapse. We also discuss in turn the impact of 
the different parameters acting upon the shell collapse and, 
thus, upon the temporal growth of the shell fragments in 
which further star formation may be stimulated. In Sect. 3, 
we show how the conditions required to stimulate a star 
formation episode within the shell provides a natural expla- 
nation to the observed metallicity range of halo GCs and 
how the shape of their metallicity spectrum constitutes the 
next step to work on. Sect. 4 describes some effects which 
our forthcoming computations should take into account in 
order to refine the present model. Finally, our conclusions 
are presented in Sect. 5. 



2 STIMULATED STAR FORMATION IN 
PROTO-GLOBULAR CLUSTER CLOUDS 

The shell will in general contain perturbed (transverse) ve- 
locity components and perturbations in the column density 
whose development leads to the transverse collapse of the 
swept PGCC and, thereby, to the formation of a second 
stellar generation. We now derive the conditions to get such 
a collapsing shell. The elementary method described below 
is adopted. 



2.1 Modelling the transverse collapse of the shell 

The computations are based on the linear perturbed equa- 
tions of continuity and motion for transverse flows in the 
shell (e.g. Elmegreen 1994). 

The perturbed equation of continuity (mass conservation) is 

where subscript T means that the gradient component 

under consideration is the transverse one, and the perturbed 
equation of motion (momentum conservation) is 



dv 
di 



-(TO 



Rs 



Cs Vcri -I- (7051 . 



(2) 



In these equations, Rs and Vs are respectively the radius 
and the velocity of the shell, ctq is the unperturbed surface 
density, ai is the perturbed surface density, v is the per- 
turbed (transverse) velocity, Cs is the velocity dispersion of 
the material inside the shell, gi is the perturbed gravity, this 
one being related to the surface density through (Elmegreen 
1994): 



gi = —2mGai 



(3) 



As mentioned above, the evolution of the shell is studied 
while propagating through the hot background, i.e. when 
the whole cloud has been swept inside the shell. The hot 
background being a diffuse medium, the shell mass Ms{t) 
does not increase any longer at this stage of its propagation 
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Figure 1. Evolution with time of the slicU radius inferred from 
Eqs. (5 - 8) assuming that 200 SNell explode at a constant rate 
during 30 million years and for three different hot protogalactic 
background pressures, from top to bottom 5 X 10~^^ (plain curve), 
10~i° (dashed-dotted curve) and SxlO"!" dyne.cm"^ (dotted 
curve). In each case, an arrow indicates the time tem at which 
the shell crosses the interface between the cold and hot phases 



and is given by the mass M of the progenitor cloud. The 
unperturbed surface density of the shell is thus given by 

1 Msit) _ J_ M 
^ 47r i?2 • 



(70 = 



47r 7?| 47r i?i ' 

Equation 1 shows that the development with time of 
any perturbation of the shell surface density (dai/dt > 0) 

is inhibited by the stretching of the perturbed region due to 
the shell expansion (i.e. Vs > 0, first term on the right hand- 
side, hereafter rhs) while the convergence of the perturbed 
flows supports the growth of the perturbation (second term 
on the rhs) . Equation 2 shows that an initial transverse flow 
of material along the shell develops (dv/dt > 0) only if the 
self-gravity (third term on the rhs) overcomes the stabilizing 
effects of the stretching (first term on the rhs) and of the 
internal pressure (second term on the rhs) , here represented 
by (?s, the shell sound speed squared. 

In order to solve Eqs.l and 2 properly, we now turn 
to the determination of the expansion law of the shell, i.e. 
Rs (t) and Vs (t) , while it propagates throughout the hot pro- 
togalactic background. 



2.2 Supershell propagation throughout the hot 

background 

The equations describing the propagation of a supernova- 
driven shell are as follow (Castor et al. 1975). 

(i) The supernova explosions add energy to the bubble 
at a constant rate Eo and the dominant energy loss of the 
bubble comes from the work against the dense shell, hence 
the variation with time of the energy Et, of the bubble obeys 

Et^Eo- 4nRs'PbRs ■ (5) 

We assume that the kinetic energy of every SN is 10^^ ergs 
and that the SN phase lasts about thirty million years, i.e. 
Eo = iV10^^ergs/30Myr. 
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Figure 2. Comparison between the evolution with time of the 
shell radius computed from Eqs. (5 - 8) and the average radius 
given by Eq. 9 (i.e. < Ra(t) >cx: t^^^) for the quoted numbers N 
of SNell and background pressures P/j 



(ii) The internal energy Ef, and the pressure Pb of the 
bubble are related through 

^R/Pt = lEb, (6) 

(iii) The shell motion obeys Newton's second law 

j^lM^mS)] = 47vR.\P, - P^^t) - . (7) 

where Pext is the pressure of the medium just outside the 
shell. 

(iv) Considering the case of swept PGCCs propagating 
through the hot protogalactic background, the mass of the 
shell is constant in time and is given by the mass of the 
cloud (soo Sect. 2.1): 



Ms(t) = M 



(8) 



The pressure external to the shell is exerted by the sur- 
rounding hot background (i.e. Vext = Ph in Eq. 7) and is 
therefore assumed to be constant in time. This pressure will 
strongly decelerate the shell as illustrated below. Numerical 
resolution of Eqs. (5 - 8) provides Rs{t) and thereby Vs{t) 
and <Jo{t). The initial conditions are those at the time tem, 
i.e. when the shell crosses the interface between the cloud 
and the surrounding background: the mass and radius of the 
shell are those of the pressure-bound cloud {Rs{tem) = R, 
M sit em) — M), the velocity of the shell is determined by its 
former propagation at constant speed V through the cloud 
(see Paper I, Eq. 13), that is Vs = and As = 0. 

Figure 1 shows the evolution with time of the shell ra- 
dius for 200 SNell, i.e. the maximum number of supcrnovae 
that a PGCC can sustain (disruption criterion, Paper I), 
and for 3 different values of the hot protogalactic background 
pressure, i.e. Pn = 5x10"", 10"" and 5x10-^° dyne.cm'^ 
The SN rate being the same in the three cases, the shells 
propagate at the same velocity in the cold phase (Eq. 13 in 
Paper I). Their expansions begin to differ once they have 
crossed the interface between their respective PGCC and 
the hot background in which the latter is embedded, this 



time being indicated by an arrow in Fig. 1. Obviously, the 
propagation of the shell through the cloud is much shorter 
than its propagation through the background. A significant 
part of the shell expansion through the background takes 
place at early time. Indeed, after a transient phase during 
which the velocity of the shell does not differ markedly from 
its velocity inside the cloud (Fig. 1), the overall expansion 
slows down and the radius of the shell scales roughly as t^^^ . 
This average expansion can be obtained from Eqs. (5 - 8) 
assuming that Ms(t) = 0, reflecting thereby that the tem- 
poral evolution of the shell radius does not depend strongly 
on the mass (Brown et al. 1995): 



1/3 



1/3 



(9) 



Figure 2 shows the good agreement between Eq. 9 
and the result of the numerical integration over time of 
Eqs. (5 - 8). Equation 9 shows that, during the long- 
term evolution, the expansion rates of two shells hav- 
ing the same [N/Ph) ratio are similar. Figure 2 illus- 
trates this effect for two sets of values, namely A'^=100 
and Ph,=5xl0~^' dyne.cm"^ (plain curve) and A''=200 and 
Pa=10~^° dyne.cm"^ (dashed-dotted curve). 



2.3 Growth with time of the shell fragments 

In order to assess whether the shell transverse collapse pro- 
ceeds successfully or not, we now numerically integrate over 
time Eqs.l and 2 in order to derive the temporal evolutions 
of the perturbed and unperturbed surface densities, d\{t) 
and cro(t), respectively. 

Assuming that the perturbed quantities follow a com- 
plex exponential of the angular position (j) along the shell, 
we get: 

-iri<l} 



(10) 



o"i(t,<A) 

and 

w(t,0) = v(t) e"''''* e*'^'*. (11) 

In these equations, Ai^ represents the phase difference 
between the perturbed surface density ai and the perturbed 
velocity v. q is the angular wavenumber and is related to the 
spatial wavenumber k by: 

r] = kRs = ^ Rs (12) 



where A is the wavelength of the perturbation, namely the 
average distance between forming fragments, the sites of fu- 
ture star formation. Therefore, r) is the number of forming 
clumps along a shell circumference and, as such, r] must be 
an integer. Moreover, any realistic perturbation must fit in- 
side a fraction of the shell circumference, say, \ < Rs or 

n > 6. 

At a time t and an angular position along the shell, 
the shell surface density era obeys 

(Ts{t, (t>) = (7o{t) + ai{t, 0) = (7o{t) + (fi{t) cos{r](t>) . (13) 

The fragmentation of the shell, i.e. its fully developed 
transverse collapse, occurs when 



cri(t) = iTo(i)- 



(14) 
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Writing —irj/Rs for the gradient and using Eqs. 3 and 
4, Eqs. 1 and 2 successively become: 

> Perturbed equation of continuity: 



ddi „ 14 . iri , 

ddl _ nVa , i n M V iA«l 

dt - ^Rs "'^ 47r e 
> Perturbed equation of motion: 



(15) 



iA,/> dv 



dv 
'di 



-ao — V e ^ +Cs (Ti 

lis 

-2mG (To <Ti 



4n IT] Cs . _iA^ 
H — -n-s <Ti e 



M 



-2mG(fie-'^'^ 



(16) 



Seven parameters intervene in the ability of the shell to 
undergo a transverse collapse: 

- the supernova number A'^ and the pressure external to 
the shell P^, which determine the shell expansion law {Rg 
depends on A'' and Ph) and mass (M depends on Ph), 

- the initial values of the perturbed quantities, ai{tem) 
and v{tem), and the associated phase difference A0, 

- the number rj of forming clumps embedded along a shell 
circumference, 

- the velocity dispersion Cs of the shell material which is 
related to the internal pressure Ps of the shell, supporting 
it against transverse collapse, through Ps = Cs^ps (ps is the 
volumic mass density of the shell material). 

In what follows, the influence of each of the parameters 
involved in the fragmentation process is investigated, the 
final aim being to check whether some reasonable sets of 
conditions can lead to the shell fragmentation. 



2.3.1 The phase difference between v and ai: A4> 

The amplitude of the perturbed surface density at a time 
t, ai{t), results from its initial value cri(iem) and from the 
transverse flows which redistribute the shell mass accumu- 
lated while the shell was propagating throughout the pro- 
genitor cloud. The fragmentation is therefore favoured if 
these transverse flows converge towards the clumps initially 
present within the shell. Such a situation corresponds to a 
phase difference A(j> ~ —\ between the two perturbed quan- 
tities (Ti (t, 4>) and v(t, (j)), that is, the transverse velocity ex- 
hibits a phase delay of a quarter of a wavelength with respect 
to the perturbed surface density (see Fig. 3). Replacing 
by this value, Eqs. 15 and 16 respectively become: 



'dt 

dv 
'di ~ 



Vs r,M V 



Rs 



4 n T] Cs 

M 



Rs cfi -I- 2irG di 



(17) 



(18) 



In what follows, the roles of the other parameters is 
studied assuming that A(j) = —tt/2, i.e. through the numer- 
ical integration of Eqs. 17 - 18. 
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Figure 3. Schematic description of the sinusoidal (i.e. com- 
plex exponential) behaviour with of the perturbed quantities 
V (lower curve) and cri (upper curve) along a fraction of the 
shell circumference at a given time assuming that the number 
of clumps r? is 15 and the phase difference Ac/) between v and cri 
is — ^ . Such a phase difference corresponds to the convergence of 
the transverse (perturbed) flows (indicated by the thick arrows) 
towards the initial clumps (filled circles). The open circles rep- 
resent the shell regions which are progressively depleted by the 
transverse flows 



2.3.2 The global parameters: N and Ph 

In order to assess whether the shell fragments or not dur- 
ing its propagation through the hot protogalactic back- 
ground, we compare in Fig. 4 the evolutions with time of the 
perturbed and unperturbed surface densities, namely ai{t) 
and cro{t), for different values of the external pressure Ph 
(5xl0"^° and 10"^° dyne.cm"^) and numbers iV of SNell 
(100 and 200). The corresponding metallicity is indicated in 
each panel. The other parameters are kept the same in every 
case: the sound speed Cs of the shell material is Ikm.s"^, 
the perturbed velocity v{tem) is Q.QlVa{tem), the number 
of forming clumps r/ is 10 and the initial perturbed surface 
density is assumed to be of order one per cent of the unper- 
turbed value, namely ai{tem) = 0.01 x ao{tem). This choice 
may seem rather arbitrary but the next paragraph will show 
that the initial amplitude of ai does not affect the results 
significantly . 

While the unperturbed surface density decreases due 

to the shell expansion at constant mass, the perturbed sur- 
face density grows at a rate which depends on the exter- 
nal pressure and on the number of SNell. The compari- 
son of the different panels in Fig. 4 shows that the de- 
velopment of di{t) with respect to (To(i) is favoured by a 
a high background pressure and a low number of SNell. 
Among the three panels presented (Ph = 10~^"dyne.cm 
and iV=100, Ph = 5 x 10 
5 X 10 " 



^° dyne.cm"^ and 7V=100, Ph = 



dyne. cm ^ and N=200), the transverse collapse 
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of the shell is successfully achieved^ when the lower number 
of supernovac is combined with the larger external pressure 
(Ph = 5 X 10"^° dyne.cm"^ and N=WO). This effect comes 
from the dependence of the shell surface density on the ex- 
ternal pressure Ph and the supernova rmmbcr A''. Indeed, 
combining the mass of a pressure-truncated gas cloud (i.e. 
M oc ) with Eq. 9, the surface density of the shell 

scales as: 

P4^)^ A^- ^''^ 

Consequently, the larger the pressure and/or the lower the 
supernova number, the larger the shell surface density and 
the larger the growth rate of the perturbation (see Eq. 1 
and Eq. 2, a larger surface density favours the terms pro- 
moting an efficient transverse collapse) . Therefore, although 
the presence of explodhig massive stars is required to stimu- 
late the formation of a second stellar generation, too large a 
number of SNell inhibits the ability of the shell to collapse 
and to form new stars^ . 

The influence of the surface density of the shell on its 
collapse is reminiscent of the star formation law on large 
scales, i.e. averaged over entire galactic discs. In that case, 
providing that the gas surface density is larger than a den- 
sity threshold, the star formation rate follows a power-law 
of the gas surface density (the so-called Schmidt law) while 
it falls sharply below (Kennicutt 1989). The threshold sur- 
face density varies from one galaxy to another but remains 
nevertheless in the range 10^°- lO^'^ cm~^. It may not be a 
coincidence that the shell surface densities obtained in the 
frame of this model are of the same order of magnitude or 
even larger (see Sect. 2.3.3). 

The dependence of the perturbation growth rate on N 
and Ph raises a point of interest, worthy of a mention here. 
It just so happens that the parameters determining the fi- 
nal inctallicity of the proto-clustcr, namely the number of 
supernovae {N, which determines the amount of metals dis- 
persed within the PGCC) and the background pressure [Ph, 
which determines the mass of the cloud, that is, the mass of 
primordial gas to be chemically enriched) are also some of 
those influencing the ability of the shell to form new stars. 
If some combinations of external pressures and SN numbers 
support the completion of the shell transverse collapse much 
more than some others do, then the formation of second stel- 
lar generations with the corresponding metallicities will be 
favoured. Therefore, in the frame of the self-enrichment sce- 
nario, there is a direct link between the achieved metallicity 
and the probability of forming halo stars, i.e. the study of 

^ At this point, it should be kept in mind that once the shell has 
achieved its complete transverse collapse, Eqs. 1 and Eqs. 2 are 
no longer valid. Therefore, once ao < cri, the curve cri(t) has no 
longer physical meaning. 

^ Unlike the protogalactic shells we study, the collapse of shells 
expanding in the Galactic disc is made easier by a larger num- 
ber of SNcII. In fact, assuming that the surrounding interstellar 
medium is roughly homogeneous, i.e. the shell radius is smaller 
than the density scale-height of the Galactic HI layer, one gets by 
mass conservation, assuming a pre-shcU density po, ctq = . 
In this case, a larger number of SNell leads to a larger radius at 
a given time and, therefore, to larger surface density and pertur- 
bation growth rate. 



the fragmentation process may shed light on the metallicity 
distribution function of Galactic halo field stars and GCs. 

Figure 4 and Eq. 19 show that the shell fragmenta- 
tion is favoured by a low number of SNe and a high pres- 
sure of the hot protogalactic background. Now, let us imag- 
ine that, on the contrary, the collapse efficiency increases 
with both decreasing number of SNell and pressure. As 
Ph and N get smaller, the number of successful transverse 
collapses increases and the newly formed stars are more 
metal-poor. Thus, the metallicity distribution function of 
the Galactic halo would exhibit an increasing number of 
proto-clusters/halo stars with decreasing metallicity. If, on 
the other hand, the shell ability to achieve fragmentation 
increased with both increasing external pressure and SN 
number, then the larger Ph and N, the more numerous 
the transverse collapses and the more metal-rich the newly 
formed stars. As a result, there would be an increasing num- 
ber of proto-clusters/halo stars with increasing metallicity. 
Neither an increasing nor a decreasing metallicity distri- 
bution function is observed for the Galactic halo. In con- 
trast, the metallicity distributions, for both halo field stars 
(Laird et al. 1988) and halo GCs (Zinn 1985), are peaked- 
shape. While the finding that the shell transverse collapse is 
favoured by large external pressures (promoting "large" , i.e. 
mildly metal-poor, metallicities) and low SN numbers (sup- 
porting low metallicities) is not sufficient in itself to draw 
some definitive conclusions regarding the shape of the halo 
metallicity distribution, it appears that it does not contra/- 
dict it either. 

2.3.3 The initial perturbed surface density: ai{tem) 

Protogalactic shells such as those studied here exhibit num- 
ber surface densities of order several 10^" cm~^. Indeed, us- 
ing Eqs. 4 and 9, one gets: 

ao{t,N,Ph) ~ 5xlO''^P^(%iV2-?/^t-2/3g.cm-2 

^ 2.5 X 10^^P,\%iV,-ort6 ^/^m- (20) 

where the subscript (10) means that the pressure is 
expressed in units of 10~^" dyne.cm"'^, the subscript 200 
means that N is expressed in units of 200 SNe and the 
subscript 6 means that the time is expressed in units of 
10® years. Considering the interstellar medium of the Galac- 
tic disc with roughly the same surface density (though a 
bit lower, i.e. ~ lO'^^cm"'^), the relative non-uniformities 
in the number surface density are typically of order 0.01 - 
0.2 (Wunsch & Palous 2001). However, the initial perturba- 
tion in the shell surface density is certainly not as large as 
some of the clumps of this quiescent disc interstellar medium 
since turbulent mixing is expected to take place within the 
shell (this is an important requirement to achieve the effi- 
cient mixing of the supernova heavy elements with the cloud 
gas and, therefore, to maintain the chemical homogeneity of 
the shell, that is, of the proto-globular cluster; see Brown et 
al. 1991), thus lowering the surface density inhomogeneities. 

Figure 4 also shows the temporal evolution of ao and 
5"! assuming that the initial perturbed surface density is 
lower than in Sect. 2.3.2, namely 0.25% a'o(iem) instead of 
l%cro(tem). Despite these different initial values, each panel 
displays almost identical ai (t) curves in both cases. Owing 
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Figure 4. Temporal evolution of cto (plain curves) and <Ji 
(dashed-dotted curves), expressed in units of IMq.pc"-^, con- 
sidering two different initial perturbed surface densities (see the 
keys) and the external pressures and SNII numbers indicated. In 
each panel, Cs=lkm.s~^, r; = 10 and v{tem) = 0.01Vs{tem.) 



2.3.4 The initial perturbed velocity: v{t.am) 

All the results presented here above assume the spherical 
symmetry of the supcrsholl. It is clear however that such a 
system is not expected to remain perfectly spherical. The 
deviations from spherical symmetry can arise, for instance, 
from the non point-like nature of the energy input, namely 
all the stars of the first generation cluster are not located 
exactly at the centre of the cloud/shell. A rough estimate 
of the initial amplitude of the transverse motions within the 
shell can therefore be derived from the size of this cluster. A 
star located at a distance r from the shell centre will induce 
a transverse velocity v such that: 



r 



(21) 



In order to estimate the size r of the cluster of massive 
stars hosted by a PGCC, we refer to R136, the dense core 
of the 30 Doradus Nebula located in the Large Magellanic 
Cloud. The 30 Doradus nebula shows an impressive exam- 
ple of a two-stage stellar formation. The energetic activity 
of a very compact bright cluster, R136, which includes sev- 
eral tens of O stars, triggers the formation of a new stellar 
generation revealed by numerous infrared sources in or near 
some bright filaments west and northcEist of R136 (Rubio 
et al. 1998). Based on Hubble Space Telescope photome- 
try, Campbell et al. (1992) detected about 160 stars more 
massive than 10 in R136 which they define as a region 
of 2.2 pc X 1.9 pc. This number of massive stars being re- 
markably similar to the numbers of SNell used in our self- 
enrichment model, we adopt R136 as the most similar ex- 
ample in the Local Group of what may have been the first 
generation cluster. A radius of 1 pc appears therefore as a 
reasonable estimate of the size of this cluster, the source of 
the energy input. 



10- 



Considering an average background pressure of 
'^'^ dyne.cm"^ and the corresponding cloud radius (~ 



30 pc, see Eq. 3, Paper I), the transverse velocity when the 
shell reaches the cloud boundary is: 



Vitem) 



Ipc 



em / - 



(22) 



30 pc 

Figure 5 displays the evolution with time of ao and ai 
assuming three different values for the initial transverse ve- 
locity, namely O.OlVsitem), 0.02Vs{tem) and Q.OSVsitem), 
while keeping all the other parameters to their previous 
values. In sharp contrast with ai{tem), v{tem) appears to 
be an important parameter of the shell transverse col- 
lapse. For instance, considering the bottom panel in Fig. 5 
(P^ = 5 X 10~^° dyne.cm"^ and iV=200), the fragmentation 
takes place less than 15 million years after the first SN ex- 
plosion if v{tem,) = 0.03Vs{tem), whcrcas it is prevented if 

V{tem) = 0.01Vs{t em I- 



to the initial transverse flows, the initial clumps quickly un- 
dergo a replenishment in shell material leading afterwards 
to very similar temporal evolutions of the perturbed surface 
density whatever the initial perturbed surface density. Ob- 
viously, this one is not a key parameter of the fragmentation 
process. 



S.3.5 The number of forming clumps: rj 

In supernova-driven shells of gas, star formation does not 
take place all along the whole periphery but, instead, takes 
place in regularly spaced clumps (i.e., where the transverse 
collapse brings some gas from depleted adjoining regions in 
the shell, see Fig. 3). The number rj of such clumps along 
a shell circumference is thus related to the wavelength A 
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Figure 5. Temporal evolution of cto (plain curves) and cti 
(dashed-dottcd curves), expressed in units of IMq.pc"-^, con- 
sidering three initial perturbed velocities (see the keys) and the 
external pressures and SNII numbers indicated. In each panel, 
Cs=lkm.s~^, r;=10 and cri(tem) = 0.01(To(tem) 
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Figure 6. Temporal evolution of ctq (plain curves) and cri (dashed 
and dotted curves), expressed in units of IMq.pc"-^, consider- 
ing the external pressures, SNII numbers and numbers of form- 
ing clumps indicated (see the keys). In each panel, Cs=lkm.s~'^, 

O-l(tem) = O.OlCTo(tem), "(tern) = 0.01Vs{tem) 



of the perturbation (Eq. 12). This aspect of shell collapse 
modelling following which stars are formed in discrete stel- 
lar subsystems along a shell periphery is supported by sev- 
eral observations of star forming shells, the so-called Sextant 
being one of the most illustrative. In the Large Magellanic 
Cloud, a set of five OB associations are located along a HI 
supershell, sustaining there about 1/6 of a complete circle 



and being thus called the Sextant (Efremov & Elmegreen 
1998). Efremov, Ehlerova & Palous (1999) showed that the 
OB associations are regularly spaced, the deprojected av- 
erage distance between two subsequent stellar groups being 
~ 37 pc. Furthermore, using numerical simulations of shells 
expanding in a mass model of the Large Magellanic Cloud, 
they concluded that the formation of these 5 OB associa- 
tions is most probably the result of a triggered star forma- 
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tion episode in the supershell created by SNell located near 
the Sextant centre. Their simulations also explain how pro- 
jection effects make the visible star forming regions only a 
fraction of a total circle. 

The influence of the number of clumps embedded within 
the shell on the transverse collapse is illustrated in Fig. 6. 
The parameter 77 is not as negligeable as the initial perturbed 
surface density. In some cases, it acts upon the fragmenta- 
tion process almost as strongly as the initial transverse ve- 
locity. It is thus worth estimating the angular wavenumbers 
which are the most favourable to the growth of an initial 
perturbation. In order to do so, wc now derive an analytical 
approximation of the instantaneous growth rate of the per- 
turbation. This is done straightforwardly by assuming that 
the perturbed quantities vary exponentially with time t, in 
analogy with the exponential growth rate found in other in- 
stability problems (e.g. the Jeans mass) even though there 
is no real exponential growth in this problem owing to the 
shell expansion and the corresponding time dependence of 
fTQ. It is thus important to keep in mind that the perturba- 
tion growth rate derived below (Eqs. 28 and 31) is not used 
to compute the temporal evolution of ui (the cases displayed 
in Fig. 6 are obtained from solving Eqs 17-18) but merely 
to derive r] estimates favouring the shell collapse. The per- 
turbed quantities arc written: 



ai{t, (f>) = ai{tem) e" e 
and 



(23) 



(24) 



where uj is the angular frequency of the perturbation. 
Following Eqs. 23 and 24, we write ui for the time derivatives 
and —iri/Rs for the transverse gradients. Therefore, Eqs. 1 
and 2 become respectively 



ujai = -2 — ai -I- cro-5-v 
its Kg 

and 

V3 , 2 irj 
its its 
using Eq. 3 in the latter. 



(25) 



(26) 



The elimination of the perturbed quantities ai and v 
between Eqs. 25 and 26 provides the dispersion equation, 
namely the relation between the angular frequency uj (i.e. 
the instantaneous growth rate) and the angular wavenumber 
r] of the perturbation: 

whose solution is given by 



= -2 E: + \2E? 



■R. 



Ri 



(28) 



Let us consider the iJrst growing mode. This one corre- 
sponds to the sequence of values of 77 which maximises the 
angular frequency co at each moment of the shell propaga^ 
tion, i.e. rj/g = ri{t) such that 

duj 



dr) 



= . 



(29) 



Equation 29 indeed corresponds to a maximum since the dis- 
criminant of Eq. 28 shows a negative curvature with 77. The 
instantaneous angular wavenumber and the instantaneous 
angular frequency associated to the first growing mode obey 
respectively 

' (30) 



ttG „ 

Vfg = —f^oRs = 



Cs 

and 



4cs2 Rs 



3Vs_ 
'2Rs 



(31) 



Taking into account the dependence of M and Rg 
(Eq. 9) on Ph and A'^, Eq. 30 shows that ?7/g depends on 
the external pressure Ph, on the SN number N and on time 

t as 



-1/3,-1/3 



(32) 



A first guess of a favourable angular wavenumber can be 

estimated from the temporal average of Eq. 30 over the time 
spent by the supershell in the hot protogalactic background: 

< "^f" >= Ai-, I ^^«(*') ' (33) 

t'em J tem 

whore At is the duration of the SN phase. Among 
the cases displayed in Fig. 6 (cs=lkm.s~^, d\{tem) = 

O.Olcro(iem), V{tem) = 0.01 (tern)) , WC 866 that the shcU 

transverse collapse is achieved if, for instance, Ph = 
lO^^'dyuc.cur^, A=100 and r;=16 (top panel) or if Ph = 
5 X 10""' dync.cnr^, iV=100 and r;=10 (middle panel). It is 
interesting to note that these values of rj reasonably match 
those given by Eq. 33, for the above mentioned cases, i.e., 

< rjfg >=12 and < rjfg >=9, respectively. As Eq. 33 has 
been derived under the assumption of an exponential growth 
rate, this agreement a posteriori justifies its validity as a 
convenient estimate of the number of forming clumps in col- 
lapsing shells. 



2.3.6 The shell sound speed: Cs 

The development of a gravitational instability within a su- 
pershell also depends on the sound speed Cs of the shell ma^ 
terial. Indeed, Cs is directly related to the thermal pressure 
Pe of the shell gas through 

Ps = = psCs . (34) 

Psmn 

In this equation, k and vriH are the Boltzmann constant and 
the hydrogen mass, respectively, while pa, Tg and p,s are 
the mass density, the temperature and the mean molecular 
weight of the shell, respectively. The larger the velocity dis- 
persion, the more the layer resists gravitational collapse. In 
an HI layer (T ~ 100 K, p, ~ 1.3 mn), the sound speed 
is ~ 0.8km.s~^, while it can be as low as ~ 0.3km.s~^ in 
an H2 layer (T ~ 20K, /i ~ 2.1 mif) (McCray & Kafatos 
1987). However, the turbulence and magnetic fields of the 
shell will increase these values. To some extent, they can be 
represented by an additional pressure term in the expression 
of Cs, i.e. 

1/2 

Cs = I h -r— + i s I (35) 
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Figure 7. Temporal evolution of ctq (plain curve) and <ti (dashed 
curves), expressed in units of lM0.pc~^, considering the sound 
speeds (see the key), external pressures and SNII numbers in- 
dicated. Other parameters are ??=10, CTi(tem) = O.OliTo(tem), 
v(tem) = O.QlVsit em } 



where Bs is the magnetic field in the shell and Tg is the 
contribution of turbulence. 

Figure 7 shows how highly sensible to Cs the fragmen- 
tation process is. Considering 100 SNell and a hot back- 
ground pressure of 5 x 10^^" dyne.cm^^, the fragmeutatiou 
takes place even with a low initial transverse velocity, i.e. 
v{tem) = O.OlVs(tem), if Cs=lkm.s~^. Increasing the lat- 
ter by 25%, the transverse collapse is very weakened and 
the fragmentation is prevented. During the last ten million 
years, the evolutions with time of ao and ai are similar, that 
is, the evolution of ai {t) is mostly driven by the dilution of 
ao(t) due to the shell expansion. 

Equation 35 illustrates the difficulty of estimating the 
sound speed of a gas. It implies the computations of its cool- 
ing history, its magnetic fields and turbulence. Moreover, the 
high sensibility of the fragmentation issue to Cs (see Fig. 7) 
shows that its value must be estimated with some accuracy. 
Such a task is well beyond the scope of the present work and, 
in what follows, we adopt Cs = Ikm.s"^, in agreement with 
many studies of supershell fragmentation (e.g. Comeron & 
Torra 1994, Ehlerova & Palous 2002). 



3 DISCUSSION 

3.1 Stars with halo GC metallicities 

The previous section has shown that the shell/swept PGCC 

may become gravitationally unstable and finally break into 
fragments providing that some conditions are fulfilled, e.g. 

V{tem)/Va{tem) — 0.03, Cs ~ Ikm.S"^, TJ ~< r]fg >. This 

is not to claim that all supershells will encounter such 
favourable circumstances, but one may expect that at least 
some of them will do. 

Figure 8 presents the results of shell fragmentation sim- 
ulations for varying values of TV, Ph, r) and v{tem)/Va{tem), 
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Figure 8. Summary of the {N , P^, v{tem)/Vs{te.m), v) values 
which may (open or filled circles), or may not (no symbol), lead 
to a swept PGCC fragmentation providing that the transverse 
flows converge towards the clumps embedded within the shell, i.e. 

= — 7r/2. The sound speed and the initial perturbed surface 
density are Cg = Ikm.s"^ and cTi(tem) = 0.01 cr()(tem), respec- 
tively. The top, intermediate and bottom panels correspond to 
initial transverse velocity of 1 , 2 and 3 per cent of the shell veloc- 
ity when it emerges out of the cloud, respectively. The isometal- 
licity curves corresponding to [Fe/H]=— 1.2, —1.5, —2, —2.5 are 
also plotted, giving thus the metallicity achieved through self- 
enrichment for each parameter combination. For eaich {N,Ph) 
pair, 54 couples of initial perturbed velocities {v{tem.)/Vs{tem)= 
0.01, 0.02, 0.03) and clump numbers (r; ranging from 6 to 40 by 
step of 2) are tested. In each panel (i.e., for each value of the ini- 
tial perturbed velocity), a couple (A'', P^) leading to a successful 
transverse collapse is marked by an open/filled circle as well as 
by the range of rj values leading to fragmentation. The larger the 
range of r] values, the bigger the symbol (see text for details) 
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assuming that Cs — l km.s ^ and ai{tan)= 0.01 CTo(iem)- Five 
hot protogalactic background pressures {Ph — 10~^^,3.2 x 
10-",10"^",3.2 X 10~^°, 10"'' dyne.cm"^) and 3 SN num- 
bers (A''=50, 100, 200) are tested. The upper limit for N is 
the maximum rmmbor of supornovao that the GC gaseous 
progenitor can sustain, namely A''=200 (i.e., if A'' > 200, 
the absolute value of the cloud binding energy is lower than 
the shell kinetic energy: disruption criterion, Paper I). The 
lower limit is imposed by c^^"-"^ , the sound speed of the 
PGCC material. Indeed, the shell is "built" while sweeping 
the PGCC and such a meiss accumulation into a shell re- 
quires the velocity of the shell to be larger than the sound 
speed of the ambiant medium. Therefore, the lower limit to 
the shell velocity in the PGCC obeys: 



yPGCC ^ ^PGCC = 8.3km.s-^ (36) 

where T and /x are the temperature (~ 10* K) and the mean 
molocular weight (~ 1.2) of the PGCC, respectively (Fall 
& Rees 1985). The number of SNell corresponding to this 
lower limit of the shell velocity is AT ~ 50 (Eq. 13, Paper I). 
Regarding the upper value of Ph, we refer to Murray & Lin 
(1992) who showed that the hot protogalactic background 
pressure depends on the galactocentric distance D as 

Pft = 1.25 X 10"^£)fcp^^ dyne cm-^ (37) 

Thus, Ph is of order 10"^ dyne.cm"^ in the very inner Galac- 
tic regions (i.e., D ~ Ikpc). 

In order to see what is the metallicity achieved in the shells 

which succeed in forming new stars, iso-motallicity curves 
corresponding to [Fc/H]=— 1.2, —1.5, —2 and —2.5 (i.e., 
metallicities typical of the Galactic halo GCs) are displayed 
in each panel of Fig. 8. 

For each couple {N,Ph), 54 combinations of v{tem)/Vs{tem) 
and ri (i.e. 3 initial perturbed velocities x 18 numbers of 
clumps) have been run. The initial transverse velocities cor- 
respond to 1, 2 and 3 per cent of the velocity of the shell 
when it enters the hot background and the number of clumps 
ranges from 6 to 40 by step of two. It appears that the shell 
is unable to fragment if r; > 40. When a set of conditions 
leads to a successful fragmentation, the corresponding point 
in the {N, Ph) diagram is marked by a circle as well as by the 
range of r] values leading to successful transverse collapses. 
Depending on whether the shell fragmentation takes place 
for to 25 per cent, 25 to 50 per cent, 50 to 75 per cent, 
or 75 to 100 per cent of the number of r] values tested, the 
corresponding triplet (N, Ph, v{tem) /Vs{tem)) is marked by 
an open circle, a small, a median-size or a large filled circle, 
respectively. Figure 8 confirms that the probability of suc- 
cessful transverse collapse fades away with decreasing exter- 
nal pressure and increasing number of SNell. As mentioned 
in Sect. 2.3.2, a low external pressure and a high number 
of SNcII loads to a larger radius for the shell, decreasing 
thereby its surface density (Eq. 19) and its ability to col- 
lapse. 

At high pressure, i.e. Ph — 10~^ dyne.cm"^, the abil- 
ity of the shell to get fragmented is limited by too large a 
number of SNell. Figure 8 shows however that a SN number 
as large as 200 does not prevent the fragmentation. Accord- 
ingly, the largest metallicity which can be achieved through 



self-enrichment is [Fe/H]~ —1.2. On the other hand, at low 
pressure, the transverse collapse is supported by a low num- 
ber of SNell. Combining a low background pressure with the 
lower limit on N, Fig. 8 shows that the lowest metallicity 
which can be achieved is [Fe/H]~ —2.8 providing that the 
relative initial transverse velocity is 3 per cent. While this 
lower limit in metallicity is a bit uncertain, as it is achieved 
in the case of the highest intial transverse velocity only, the 
examination of the three panels in Fig. 8 clearly shows that 
a metallicity of ~ —2.5 is actually achievable. These extreme 
values ([Fe/H]~ —1.2 and —2.5) match nicely the metallic- 
ities exhibited by the most metal-rich and the most metal- 
poor Galactic halo GCs, respectively. 

At this stage, it is worth keeping in mind that the ques- 
tion addressed here above concerns the ability of the shell 
to form stars and not yet the ability of these stars to evolve 
into a stellar cluster. We will address the ability of the newly 
formed stars to form a bound cluster in a forthcoming pa- 
per and we emphasize here that, among the cases of suc- 
cessful fragmentation displayed by Fig. 8, some of the newly 
formed stars may not be able to form a bound cluster. In 
other words, the collapsed shells may be a source for both 
halo GCs and halo field stars. 

Section 2.3.2 has discussed how the self-enrichment 
model for GC formation can shed light on the metallicity 
distribution of the Galactic halo. Because the parameters 
which determine the metallicity, i.e. N and Ph, also control 
the shell surface density, there is a direct link between the 
metallicity and the probability of star formation^ . Further- 
more, the initial radius and velocity of the second generation 
stars being the radius and the velocity of the shell at the 
time of their formation, their binding depends on the shell 
expansion law and, therefore, again, on A'^ and P/j. As a con- 
sequence, the probability of getting a bound cluster from a 
shell of newly formed stars is also related to their metallicity. 
If we assume that a significant fraction of the field content 
comes from "failed GCs" (i.e. shells of stars which did not 
succeed in forming GCs), then these relations between the 
metallicity on the one hand and the probabilities of forming 
a second stellar generation and a bound cluster on the other 
hand open the way to the computation of the metaUicity 
distributions of both halo held stars and halo GCs. 

Since the formation of a GC requires one more condi- 
tion than the shell transverse collapse (i.e. the binding of 
the second stellar generation), only a subset of the {N,Ph) 
combinations giving rise to shell fragmentation may lead to 
the formation of a bound cluster. Such a conclusion does 
not contradict, and even fits, the fact that the metallicity 
spectrum of halo field stars is larger than the one of GCs 
(e.g. Laird et al. 1988). 



^ This is not to say that the probability of star formation is deter- 
mined by the metallicity, as it can be assumed in pre-enrichment 
models in which the metallicity may control the star formation 
through line cooling processes. It rather means that the metal- 
licity and the probability of star formation derive from the same 
parameters, namely the number of SNell and the pressure of the 
hot protogalactic background 
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3.2 The very metal-poor stars ([Fe/H] < -2.5) 

The most striking difference between the halo GC and halo 
field metallicity distributions resides in the much more ex- 
tended metal-poor tail of the field compared to the GC sys- 
tem. In fact, while the most metal-poor GCs show [Fe/H]~ 
—2.5, the Galactic halo hosts field stars whoso metallicity is 
much lower (e.g., CS22876-032, [Fe/H]~ -3.7, Norris, Boors 
& Ryan 2000). The recent discovery of HE0107-5240, the 
most metal-deficient star ever discovered ([Fe/H]~ —5.3, 
Christlieb et al. 2002) has decreased even more the lower 
limit of the observed rnotallicity distribution of halo field 
stars. Our model does not seem to be able to explain these 
very metal-poor stars (see Fig. 8) and we must therefore 
investigate alternative scenarios for the formation of these 
stars more metal-deficient than the most metal-poor GCs. 

3.2.1 Isolated Type II Supemovae 

The existence of very metal-poor stars has often been at- 
tributed to star formation episodes triggered by individual 
SNell exploding in primordial gas clouds (e.g., Shigeyama 
& Tsujimoto 1998, Argast ot al. 2000, Karlsson & Gustafs- 
son 2001), that is, a formation scenario similar to the self- 
enrichment scenario for GCs except regarding the considered 
number of massive stars. The metallicity of such stars is de- 
termined by the ratio of the mass of metals ejected by the 
SNII to the mass of hydrogen gathered by the shock wave. 
The SNII yields computed by Woosley & Weaver (1995) for 
zero metallicity stars show that a SNII with progenitor mass 
m ejects a mass of metals such that ~ 0.3m — 3.5 (in 
unit of one solar mass). On the other hand, the mass Msw 
of interstellar gas collected by the shock wave is related to 
the explosion energy Eq through 

M..=5.1xl0^Me(^^) , (38) 

assuming a sound speed of lOkm.s"^ (or T ~ 10* K) for 
the interstellar gas (Shigeyama & Tsujimoto 1998). Using 
these relations and assuming Eo = lO''^ ergs for a canoni- 
cal SNII, the metallicity achieved in shells of gas driven by 
isolated SNell is straightforwardly derived (see the second 
column of Table 1). The explosion energy being the same 
whatever the SNII progenitor mass, every SNII sweeps the 
same amount of primordial circumstellar gas, irrespective of 
the SNII mass (Eq. 38). Therefore, the final metallicity is 
determined solely by the mass of heavy elements released 
in the interstellar medium by the exploding star and, thus, 
by the SNII progenitor mass. As a consequence, the more 
massive the SNII, the larger the final metallicity. As already 
noticed by Shigeyama & Tsujimoto (1998), the formation of 
very metal-poor stars with [Fe/H] ~ —4 and [Fe/H] ~ —2 
may therefore be ascribed to "low"-mass (i.e., M ~ 12 Mq) 
and high-mass (i.e., M ~ 40 M©) SNII, respectively. In or- 
der to dig out the potential link between very metal-poor 
stars and isolated SNell, Shigeyama & Tsujimoto (1998) 
also exploited the yields of core-collapse SN predicted by 
nucleosynthesis calculations (e.g., Woosley & Weaver 1995, 
Tsujimoto et al. 1995) . They focused on the abundance pat- 
terns of C, Mg, Si and Ca. The yields of these elements 
being not (or at least less) affected by the mass-cut issue, 
they should therefore be known with a better accuracy than 



the yields of heavier elements such as the iron peak ones. 
They noticed the good agreement between the abundance 
patterns arising in the remnants of first generation SNell, 
as predicted by nucleosynthesis calculations, and those in- 
ferred from the spectra of halo stars with [Fe/H] < —2.5. 
As a result, Shigeyama & Tsujimoto (1998) suggested that 
very metal-poor stars actually formed out of interstellar gas 
which had been swept by isolated SNell. They also note 
that the large observational scatter in the abundance ratios 
among these stars can be explained by the differences in 
SNII yields with the mass of the progenitor. 

It is worthy of a mention that such isolated SNell are 
certainly not able to trigger the formation of GCs owing 
to the small amount of interstellar gas swept by the blast 
wave (Eq. 38; also, not all the gas will be converted into 
stars). Should such low-mass halo clusters have managed to 
form however, they will have been quickly destroyed by the 
Galactic tidal fields, low-mass clusters being among the most 
vulnerable in this respect (e.g., Gnedin & Ostriker 1997). 
Therefore, this extension of the self-enrichement scenario to 
single SNII explosions does not contradict the absence of 
Galactic halo GCs with [Fe/H] < -2.5. 

3.2.2 Hypemovae 

In fact, even stars as massive as 40 Mg may trigger the for- 
mation of stars with a metallicity as low as —4, providing 
that they explode as hypemovae, i.e., supernovae character- 
ized by explosion energies of order Eo ~ 10®^ — 10^^ ergs. 
This class of objects includes two categories depending on 
the progenitor mass, namely core-collapse hypcrnovac and 
pair-instability hypemovae. They both have been invoked 
as possible explanations to peculiar abundance patterns ob- 
served in some very metal-poor field stars. 

Core-collapse hypemovae. Recent observations suggest 
that at least some core-collapse supernovae explode with 
explosion energies ten to one hundred times higher than 
the energy released by a canonical SNII (e.g., Galama et 
al. 1998). They likely originate from relatively massive star 
(M> 25 Mq). In contrast to pair- instability hypemovae (see 
below) , the mass range of their progenitor is thus similar to 
the one of canonical SNell. 

If hypemovae occured in the early stage of the Galactic 
evolution and induced star formation, their abundance pat- 
tern may still be observable in the atmospheres of some low- 
mass halo stars. As for the case of star formation triggered 
by single SNell, we can derive an estimate of the metallicity 
of these halo stars using Eq. 38 and hypernova yields. Naka- 
mura et al. (2001) have investigated in detail the nucleosyn- 
thesis of core-collapse hypemovae and compared it with the 
yields of canonical (i.e., Eq ^ lO^'^ ergs) core-collapse SNell 
with similar progenitor mass. Their Tables 2-5 show that 
SNell and hypemovae with similar progenitor mass release 
much the same amount of metals. The resulting metallicity 
of stars formed in hypernova remnants are given in Table 
1, following the same of line of reasoning as in the previous 
section. Hypemovae being much more energetic than SNell, 
and the mass of interstellar gas swept by the blast wave be- 
ing roughly proportional to the explosion energy (Eq. 38), 
they will collect a larger amount of interstellar gas. There- 
fore, should core-collapse hypemovae be able to trigger the 
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Table 1. Dependence of the metallicity [Fe/H] on the progeni- 
tor mass m and explosion energy Eq of supernovae (SN; Eq = 
10^1 ergs) and hypernovae (HN; Eq = 10 or 100 x lO^^ergs). 
Every mass is expressed in units of one solar mass. (1) A 12 Mq 
star explodes as a canonical SNII. According to Nakamura ct 
al. (2001), stars less massive than 25 Mq do not explode as core- 
collapse hypernovae and such cases are thus not considered in this 
Table. (2-3) More massive stars explode either as core-collapse 
supernovae or core-collapse hypernovae. (4) Very massive objects 
(i.e., M > 100 M0) explode as pair-instability hypernovae. Msw, 
the mass of interstellar gas collected by the explosion blast wave, 
depends on the explosion energy through Eq. 38. m^ is the mass 
of metals released by an exploding massive star with progenitor 
mass m. The metallicity [Fc/H] is derived assuming that the mass 
niz of metals is mixed with the mass Msw of collected interstellar 
gas 



Eo/lO^i ergs 


1 (SN) 


10 (HN) 


100 (HN) 




5.1 X 10* 


5.1 X 10^ 


5.1 X 10^ 


(1) (m,m^)=(12,0.1) 
[Fe/H] ~ 


-4 






(2) (m,m^)=(25,4) 
[Fe/H] ^ 


-2.4 


-3.4 


-4.4 


(3) (m,m2)=(40,10) 

[Fe/H] ~ 


-2.1 


-3.1 


-4.1 


(4) (m,m^)=(200,100) 
[Fe/H] ^ 




-2.0 


-3.0 



formation of new stars, these stars will be more metal-poor 
than the ones formed in the remnants of canonical SNell. 
Thus, they may be well-suited to explain the formation of 
stars with a metaUicity as low as [Fc/H] ~ —4 (see Table 1). 

The most significant feature of hypcrnova nucleosyn- 
thesis is their iron production, this one being larger than in 
SNell by a factor 2 to 10 (Nakamura et al. 2001). This leads 
to small abundance ratios of a elements over iron. Naka- 
mura et al. (2001) thus suggested that the gas out of which 
the very metal-poor binary CS 22873-139 ([Fo/H] = -3.4) 
formed had been contaminated by the ejecta of an hypernova 
as this halo star shows almost solax [Mg/Fe] and [Ca/Fe] rar 
tios. On the other hand, while stars with —2.5 < [Fe/H] < 
show [Zn/Fe] ~ (e.g., Primas et al. 2000), this abundance 
ratio is steadily increasing towards [Zn/Fo] ^ 0.5 as the 
metallicity decreases, for stars more metal-poor than [Fe/H] 
~ -2.5 (e.g., Cayrel et al. 2003). Umeda & Nomoto (2002) 
notice that such a large [Zn/Fe] ratio arises naturally in their 
own hypernova model and thus conclude that core-collapse 
hypernovae are likely to have contributed to the early Galac- 
tic chemical evolution. 

Pair instability hypernovae. These very massive (~ 
140-260 Mq) objects exploding with E > lO'^^ ergs, they 
are also called hypernovae. Their name (i.e., "pair instabil- 
ity") refers to the electron-positron pair instability process 
encountered during the central oxygen-burning stages (see 
Umeda & Nomoto 2002, their appendix for a detailed de- 
scription). 

The main feature of these very massive stars is their 
ability to "pass carbon and oxygen from the helium-burning 



core through the hydrogen-burning shell, in such a way that 
it is CNO processed to nitrogen before entering the hydro- 
gen envelope" (Carr, Bond & Arnett 1984). These stars thus 
produce large supersolar values of N/Fe, an effect not pre- 
dicted by models of Galactic chemical enrichment based on 
stars less massive than 100 Mq . Following their discovery of 
CS 22949-037, a very metal-poor star ([Fc/H] = -3.8) show- 
ing extreme nitrogen enhancement ([N/Fe] = 2.3), Norris et 
al. (2002) suggested that this star may have been formed 
out of gas polluted and compressed by such a very massive 
hypernova. 

On the other hand, in marked contrast with core- 

coUapsc hypcrnova yields, Umeda & Nomoto (2002) quote 
that pair-instability hypernovae arc unlikely to produce 
[Zn/Fe] ratios as large as in very metal-poor stars where 
[Zn/Fe] ranges from solax up to ~0.5 (Cayrel et al. 2003). 
Actually, the abundance ratio derived from their hyper- 
nova yields is of order [Zn/Fe] = — 1.5, that is, 1 to 2 dex 
smaller than in very metal-poor stars. Additionally, Heger 
& Woosley (2002) note the absence of the r-process in these 
massive objects, which conflicts with observations showing 
appreciable amounts of r-process elements in very metal- 
poor stars (e.g., Burris et al. 2000). Based on these two 
arguments, i.e., the low abundance ratio [Zn/Fe] as well 
as the absence of r-process elements in pair instability hy- 
pernova ejecta, Umeda & Nomoto (2002) and Heger & 
Woosley (2002) conclude that the abundances of very metal- 
poor stars cannot be ascribed to pair instability hypernovae 
only and must include the contribution of an additional 
nucleosynthetic component, namely those lower mass stars 
that make "regular" supernovae. We also note that pair- 
instability hypernovae release fairly large amounts of met- 
als in the interstellar medium, of order 100 Mq (Umeda & 
Nomoto 2002, their Tables 15-18). In fact, these stars dis- 
rupt completely when exploding, leaving no compact rem- 
nant after the explosion (i.e., no issue of "mass-cut" or 
"fall-back", see Heger & Woosley 2002). The mixing of so 
large an amount of metals with a mass of primordial gas 
of ~ 5 X IO^Mq (Eo = 10 X 10^^ ergs) or ~ 5 x lO'^M© 
(Eo = 100 X 10®^ ergs) will lead to metallicities of order 
[Fc/H] ~ -2 or [Fo/H] ~ -3, respectively. Therefore, the 
metallicity of the most metal-poor stars in the Galactic halo 
cannot be explained by pair-instability hypernovae. 



3.2.3 Single short-acting source of energy and triggered 
star formation 

The issue of whether isolated exploding massive or very mas- 
sive stars can stimulate the formation of now stars in the 
layers of interstellar gas they have swept is not explicitely 
addressed in the papers mentioned in Sect. 3.2.1 and 3.2.2 
(e.g., Shigeyama & Tsujimoto 1998, Argast et al. 2000, Nor- 
ris ct al. 2002). It is merely assumed that such stimulated 
star formation actually took place during the early Galactic 
stages. 

The model presented in this paper (see sect. 2) is not 
the most appropriate to address this issue quantitatively as 
it deals with a continuous input energy and not with a single 
short-acting source of energy. Actually, supershells created 
in connection to single explosion and those created around a 
cluster of massive stars do not show the same expansion rate 
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with time (see below). We have seen in Sect. 2 how the ra- 
dius and the velocity of the shell affect its ability to collapse 
through the shell surface density and the stretching of the 
forming perturbation, respectively. Hence, to solve the issue 
of whether a single supcrnova/hypernova triggers the forma- 
tion of new stars requires to derive a new temporal evolu- 
tion of the shell radius, in the PGCC and in the protogalac- 
tic background, adequate to single explosion. This is beyond 
the scope of the present paper whose main goal is to address 
the formation of stars with —2.5 < [Fo/H] < —1 (i.e., the 
metallicity range of halo GCs). Some qualitative inferences 
can however be made. Efremov, Ehlerova & Palous (1999) 
compared the propagation in an homogeneous medium of 
both kinds of supcrshclls. In case of an abrupt energy input, 
the early expansion of the supershell proceeds at a higher 
velocity than in the case of a continuous input of the same 
energy amount. Accordingly, in the former case, the shell of 
swept gas is more strongly stretched and this hampers the 
formation of an initial perturbation. After a few millions 
years however, the shell created in connection to an abrupt 
energy input slows down and shows both a radius and a 
velocity lower than in case of a continuous energy input. 
In contrast to the initial shell expansion, this second stage 
favours the transverse collapse through a larger shell sur- 
face density and a weaker stretching of the shell perturbed 
regions. Whether the collapse will occur is actually not cer- 
tain. The source of the energy input being point-like (i.e., 
its size is much smaller than the one of a cluster of mas- 
sive stars), the initial amplitude of the transverse motions 
within the shell will be lowered accordingly, thus reducing 
the collapse. 

Comparing SNell and hypernovae, hypernovae sweep 
a larger amount of gas thanks to their larger energy in- 
put (Eq. 38). However, this also leads to a larger shell 
radius. Both effects acting on the shell surface density in 
opposite ways, only detail computations will toll us which 
class of objects is the most efficient in forming new stars. A 
So ~ 2 X 10®^ ergs hypernova being able to collect amount of 
gas as large as 10® M© (Eq. 38), one might think that pair- 
instability hypernovae could be related to the formation of 
GCs, at least the most metal-poor ones as such hypernovae 
can chemically enrich the initially pristine gas up to [Fe/H] 
~ -2 (see Table 1). However, the [Zn/Fe] abundance ratio 
predicted for hypernova ejecta (i.e., [Zn/Fe] < —1.5, Umeda 
&: Nomoto 2002) shows a sharp discrepancy with the roughly 
solar value observed in halo field stars with [Fe/H] ~ —2 
(Primas et al. 2000) as well as in NGC6397, a halo GC with 
the same metallicity (Thevenin et al. 2001). Therefore, the 
abundance patterns of stars in this metallicity range cannot 
be ascribed to pair-instability hypernovae. 

Obviously, it is not straightforward to conclude whether 
single massive star explosions are able to trigger the forma- 
tion of new stars. More computations are required to answer 
a question which is certainly worthy of further investigation 
owing to its potential link with very metal-poor stars. 

3.2.4 External pollution onto PopIII stars 

Shigoyama, Tsujimoto & Yoshii (2003) have recently sug- 
gested that some very metal-deficient stars may have been 
born as metal-free (PopIII) stars whose external layers were 
afterwards polluted by the accretion of chemically enriched 



interstellar gas. The metallicity achieved by the external lay- 
ers will depend on the metallicity and mass of the accreted 
gas, as well as on the stellar mass fraction with which the 
accreted gas has been mixed. As a consequence, the chemi- 
cal enrichment of stellar superficial layers will show up more 
markedly in main sequence stars than in red giants, the con- 
vective envelope (i.e. the mixing zone) being almost two or- 
ders of magnitude larger in halo giants (about half the star 
mass) than in non-evolved metal-poor stars (about one per 
cent of the star mass). Sliigcyama et al. (2003) have shown 
that the extremely low metallicity ([Fe/H]~ —5.3) of the gi- 
ant HE0107-5240 (Christheb et al. 2002) may originate from 
such a mechanism, that is, an external pollution while the 
star was on the main sequence followed by the dilution of 
the accreted material as the star started ascending the red 
giant branch. Owing to their thin convective envelopes, pol- 
luted Pop HI stars still on the main sequence would exhibit 
larger metallicities, i.e. of order [Fe/H]~ —3, although still 
lower than the most metal-poor GCs. 



3.2.5 Star formation in dwarf galaxies 

Within the frame of a hierarchical model for halo forma- 
tion. Cote et al. (1999) proposed that the excess of very 
metal-poor field stars with respect to GCs was formed 
in the most metal-deficient dwaxf galaxies trapped within 
the Galactic potential well. These galaxies are not ex- 
pected to contribute to the GC population. Indeed, due to 
the luminosity-metallicity correlation observed among dwarf 
galaxies (Gilmore 2000, Mateo 2000), the most metal-poor 
dwarfs are also expected to be the least luminous. Dwarf 
galaxies exhibit a luminosity threshold for hosting a GC 
system. For instance, the Fornax dwarf spheroidal galaxy 
is the faintest (M„ ~ —12.3, Harris 1991) galaxy known to 
harbour its own GC system. Therefore, the dimmest dwarf 
galaxies, which are also the most metal-poor, may have con- 
tributed to the halo field and not to the halo GC system. 



4 CAVEATS AND FUTURE WORK 

It is worthy of a mention that the model results presented 
above constitute a first and preliminary step in our com- 
putations of globular cluster star formation. In fact, several 
aspects are still to be worked out. We now give some ex- 
amples of these, which our future developments will have to 
encompass. 

- As the shell collapse is followed until the perturbed and 
unperturbed surface densities reach similar values, a non- 
linear model would better describe the growth with time of 
the shell fragments. In this respect, we note that Wunsch & 
Palous (2001) studied the transverse collapse of supershells 
propagating in the Galactic disc using both linear and non- 
linear models. The comparison of the results obtained in 
both cases, and for instance the time at which the shell gets 
fragmented, does not show significantly different results. 

- At some moments of its propagation through the hot 
protogalactic background, the shell undergoes transient ac- 
celerations (see Figs. 1 and 2) which may lead to Rayleigh- 
Taylor instabilities. Rayleigh- Taylor instabilities take place 
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when a cold, dense gas (i.e. the supershell) is accelerated 
by pressure from a hot more rarefied gas (i.e. the bubble). 
Spikes and globules of the colder fluid will tend to penetrate 
the warmer fluid (Spitzer 1978). While this effect may dis- 
rupt the shell and, thus, prevent the process of triggered star 
formation, it can also lead to local overdensities in the shell 
material, accelerating thereby star formation. It will there- 
fore be of interest to estimate whether the shell transient 
accelerations may trigger Raylcigh- Taylor instabilities, how 
far these ones develop and how the process of stimulated 
star formation within the shell is affected. 

- The computation of the shell expansion law can still 
be improved as it does not yet include energy losses, such 
as those induced by viscous processes. For given number 
of SNell and external pressure, these effects may lower the 
radius of the shell with respect to what has been com- 
puted in Sect. 2.2 and, thus, increase the shell surface den- 
sity and the probability of star formation. If the forma- 
tion of stars may be triggered even at pressures lower than 
Ph = 10^^^ dyne.cm"^, metallicities lower than [Fe/H]~ 
—2.8 may be achieved. 



5 SUMMARY 

This paper has presented the results of simulations dedicated 
to the transverse collapse of shells of gas resulting from the 
sweeping of gaseous GC progenitors by SNoII. In these sim- 
ulations, the growth of an initial perturbation in the shell 
surface density is followed by solving the linear perturbed 
equations of continuity and motion for transverse flows in 
a spherical shell of gas. Such a collapse depends on several 
parameters, namely the number of SNcII, the background 
presure, the sound speed of the shell gas and the initial con- 
ditions of the perturbation (i.e. the number of clumps, the 
initial perturbed surface density and velocity and the corre- 
sponding phase difference). All these parameters have been 
discussed in turn. The results show that the pressure Ph of 
the hot protogalactic background {Ph ^ 10"^° dyne.cm"^, 
FaU & Rees 1985, Murray & Lin 1992) and the numbers N of 
SNell allowed by the disruption criterion (i.e. smaller than 
200, Paper I) can indeed lead to a successful shell trans- 
verse collapse, and thereby to the formation of new stars, 
assuming some reasonable initial conditions for the pertur- 
bation (see Fig. 8). The metallicities achieved in the shells 
able to collapse agrees with the metallicity range of Galactic 
halo GCs, namely, —2.5 < [Fe/H] < —1. Furthermore, while 
N and Ph determine the metallicity achieved through self- 
enrichment, they also control the probability of triggered 
star formation and the ability of these second generation 
stars to form a bound GC. Such a property is the most in- 
teresting since it opens the way to the understanding of the 
halo metallicity distribution functions, for both stars and 
clusters. 

This paper has been typeset from a T^JC/ file prepared 

by the author. 
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